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Abstract 

We study an opinion formation model that takes into account that individuals have 
diverse preferences when forming their opinion regarding a particular issue. We show that 
the system exhibits a phenomenon called "diversity-induced resonance" [Tessone et al. Phys. 
Rev. Lett. 97, 194101 (2006)], by which an external influence (for example advertising, or 
fashion trends) is better followed by populations having the right degree of diversity in 
their preferences, rather than others where the individuals are identical or have too different 
preferences. We support our findings by numerical simulations of the model and a mean-field 
type analytical theory. 

1 Introduction 

Resonance in forced dynamical systems is a topic of widespread interest with many applica- 
tions. A resonance is a maximum in a suitably defined response of the system to an external 
forcing. It usually requires the tuning of some system parameters to the optimal value. The 
simplest example is that of a forced linear oscillator whose amplitude of oscillations (the 
response) reaches a maximum when the natural frequency of the oscillator matches that of 
the external forcing. In coupled nonlinear systems, more complex resonances can appear. For 
example, they may mode-lock into states where the ratio between the individual frequencies 
are rational numbers in wide parameter regions, a phenomenon known as Arnold tongues 

It was shown in the early 80's that a resonance can also appear as a function of the 
intensity of the fluctuations, of either internal or external origin. The basic mechanism leading 
to this stochastic resonance E] is rather generic and in its simplest form requires only a 
bistable system, a sub-threshold periodic forcing and a fluctuating additive term (noise) in 
the dynamics. The forcing induces a periodic lowering of the barrier separating the two stable 
fixed points, so helping the fluctuation dynamics to overcome the barrier in one direction or 
the other. The matching for resonance occurs when half the period of the external forcing 
equals the Kramers' time. The surprising result that fluctuations can enhance the response of 
a dynamical system to external forcing has become a new accepted paradigm and there have 
been many extensions and applications pTJ, including neural systems [5], non- linear electronic 
devices j6j, sensory systems [7], social dynamics [8, 9J, etc. 

Although most of the work in this field has considered simple, low-dimensional systems, 
more recent work analyses the role of fluctuations in the response of an extended system 
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\10\ 111). A usual modelling is that of many interacting units located on the sites of a lattice, 
such that the individual responses to the forcing are modified by the mutual interactions. A 
typical assumption in this case is that all the units are identical in the sense that they all 
possess the same values for all constituent parameters and that there is some regularity in 
the network of interactions. For most applications, mostly in the biological or social sciences, 
those assumptions are not correct since some sort of diversity or variability will ineluctably 
be present. We have shown in a recent work [12] the existence of a new type of resonance as a 
function of the diversity of the system. That work focuses on bistable and excitable systems 
in which the diversity is modelled by quenched noise or, more specifically, by a parameter 
that adopts a different value for each of the units. Related work[13j has shown that diversity 
can also induce synchronised spiking in an extended, unforced, excitable system. 

Surprising at first, the fact that the right amount of diversity can enhance the response 
to an external forcing is not against our intuition. Think, for example, of a society which is 
very homogeneous in that all members of the population work on a particular economical 
field. If the economy tilts and that particular field becomes of less importance, it will have a 
big negative impact in the overall wealth of the population since individuals will not be able 
to follow the change. However, if there is some degree of heterogeneity and fractions of the 
populations work on different fields, there will be always a section that can adapt easily to 
the changing economy. The final ingredient that allows the whole society to follow the change 
is some degree of interaction by which the benefited agents can pull the others towards the 
new field. 

These ideas were put forward in ref. pj] where we presented a mathematical model that 
displays this effect of diversity-induced resonance. We considered an ensemble of globally 
coupled N bistable units Xi(t) whose dynamics is given by: 



The parameters a* are independently drawn from a Gaussian distribution of zero mean and 
variance a 2 . The value of a can be considered as a measure of the diversity. If a — 0, all 
systems are identical, whereas increasing values for a indicate a larger degree of heterogeneity. 
We consider that the external periodic forcing of amplitude A is sub-threshold for those 
systems with a, close to zero. This means that for a = the system as a whole is unable 
to display a wide response to the forcing and the collective variable X(t) = j* x i(t) 
oscillates around one of the equilibrium values X(t) = ±1 with an oscillation amplitude 
proportional to A. Imagine that the oscillation point is X(t) = +1. As a increases there 
will be a fraction of units (those with a sufficiently large, negative, value for at) for which 
the weak forcing is now sufficient to take them to the minimum X4 = — 1 . If the intensity of 
the coupling C is sufficiently large the whole system will be pulled by those units and taken 
to that minimum. Hence, the collective variable will have performed a large excursion from 
X(t) = +1 to X(t) = —1. The opposite jump from X(t) = —1 to X(t) = +1 is induced when 
the sign of the forcing is reversed as induced by those units that have a large, positive, value 
for di. 
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It is important to realise, as explained in the theoretical treatment of [14], that this 
resonance mechanism relies of the individual units having different dynamical response to 
the external forcing but that the origin of the heterogeneity in the response is not important. 
Different sources of disorder such as noise, diversity, non-regular network of connectivities, 
inhibitory couplings, etc. can be the origin of the resonance. The combined effect of noise 
and diversity has been analysed in [T51 [16] . This result has already been shown of interest in 
many fields, from complex networks [17] to cellular signalling [18]. A linear model in which 
a full analytical calculation is possible has been recently studied[T9]. 

In this paper, we show a rather different example of diversity-induced resonance in an 
opinion formation model. The model is a simple majority model with the addition of pref- 
erences in the individual choices. Those preferences vary from individual to individual and 
are the source of diversity. The interest of the paper is twofold. First, by giving an example 
which is very far away from the dynamical system described above, we want to emphasise the 
generality of the mechanism leading to the resonance. Second, we believe that the example 
has interest on its own in the field of social sciences, since it shows that an external forcing 
(imitating the effect of advertising) has a larger impact on a heterogeneous society than on 
a completely homogeneous one. This effect might be relevant when explaining the changes 
in opinion (e.g. in poll's results) motivated by an apparently small change in the external 
environment. 

The outline of the paper is as follows: in section [2] we define the model for opinion forma- 
tion and highlight its formal similarities to other well known models of spin-glass systems, 
while stressing the ingredients that, according to the general discussion, might lead to a 
resonance effect. In section [3] we present the results of numerical simulations that show the 
existence of the resonance as a function of a parameter measuring the diversity in the individ- 
ual preferences. In section |4j we introduce a mean-field theory that focuses on the collective 
variable and from which a global mechanism for the resonance can be extracted. Finally, in 
section [5] we end with brief conclusions and outlooks. 

2 Model studied 

Although the focus is very different, the model we have introduced bears many similarities 
with the random field Ising model [20]. This model has attracted much attention because of 
its interest for modelling disordered magnetic materials [21J and also spin-glasses [22]. More 
recent work has focused on the hysteresis behaviour when subjected to a slowly varying 
magnetic field [23J or the question on how the system can reach the global energy minimum 
[23] . The use of techniques of statistical physics to model social behaviour has a long tradition 
(see, e.g., [9]) and suitable modifications of the random field Ising model have been used 
already to model the dynamics of social systems. To the best of our knowledge, Galam |25| 
was the first one to model individual preferences by a random field. Stochastic resonance 
induced by fluctuating terms in the dynamics was described by Kuperman and Zanette [8]. 
Michard and Bouchaud |26j used the random field Ising model in an external field and studied 
the emergence of collective opinion shifts in a diverse population. The imitation mechanism, 
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included in this work, has also been presented as a key feature in price formation dynamics 
by Zhou and Sornette [27] . 

We consider a model for opinion formation in which the opinion on a particular topic is 
considered to be a binary variable (against or in favour of such a topic) . There are N indi- 
viduals, each one having an opinion Hi{t) = ±1, i = 1,...,N, at time t. The opinion Hi(t) of 
individual i can change due to (i) the interaction with the ki individuals in its neighbourhood 
n(i), modelled by a majority rule, and (ii) the influence of advertising, modelled as the effect 
of some external time- varying agent. Each individual has a tendency to favour one of the 
two opinions, +1 or —1. We introduce diversity in the fact that this preference for one of the 
two opinions is stronger in some individuals than in others. We model the effect of individ- 
ual preferences by a set of independent parameters 9i. They are drawn from a probability 
distribution g{9), which satisfies (9i) = 0, (OiOj) — Sij a 2 . Their influence becomes apparent 
when we spell out the evolution rules of this model: 

(i) Select randomly one individual i. Its opinion at time t is modified as : 



In words, individual i adopts the average opinion in its neighbourhood when this average 
opinion overcomes its preference 0j. This is a mechanism of social pressure weighted 
against individual preferences. For instance if 9i = 0.3 (resp. 6i = —0.3) it is necessary 
that the proportion of neighbours supporting the —1 (resp. +1) opinion exceeds 70% 
in order for individual i to adopt the majority opinion. Note that when \8i\ > 1 the 
individual will keep its preferred opinion no matter what the social pressure is. 

(ii) With probability A\ sin(il£) | , the opinion is set to 



This represents the effect of a time dependent external global field (advertising, for 
example). A is a measure of the strength of the field. This field has to overcome the 
preference 0i (weighted by a scale factor a) in order for the individual to adopt the 
value favoured by the field, i.e.: the advertising has more effect on those individuals 
whose preference coincides with the sign of the advertising. 

After these two steps have been taken, time increases by t —t t + dt = t + l/N and a 
new individual is selected again at random. The process is repeated for many cycles of the 
external forcing. We are interested in quantifying how well the system globally responds to 
the external forcing. To this end, we focus on the time evolution of the average opinion: 



In general, m(t) oscillates in time with the frequency f2 of the forcing. The amplitude of 
the oscillations of m(t) is a measure of the response to the forcing. An equivalent measure, 




Hi(t + dt) = sign [sin(Sli) + a6i] . 



(3) 




(4) 
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but more useful from the computational point of view, is the so-called spectral amplification 
factor R defined as[4]: 

R = 4A- 2 \(m(t)c- mt )\ 2 , (5) 

where (. . .) denotes a time average. 

The main result of this paper, as shown in the next sections, is the existence of a value 
of the diversity a for which the response R takes a maximum value. This resonance effect 
appears for weak forcing (A sufficiently small) and implies that the advertising has an optimal 
effect on the population when there is some degree of diversity in the preferences. The 
"microscopic" mechanism is easy to understand. In a diverse society, there is always a fraction 
of the population which is receptive to follow the external field. This fraction initiates a 
change in the opinion and then, by the imitation mechanism, the change is spread towards a 
larger fraction. If the population is not too diverse, the fraction that can follow the external 
signal is small and it is not enough to initiate the global change. If the population is too 
diverse, however, the imitation mechanism is not effective. In section [4] we will present a 
mean-field theory that offers an alternative, "macroscopic", explanation to this microscopic 
mechanism. It will be clear after the theoretical treatment that the origin of the resonance 
can be traced to the lack of order caused by the diversity and, as stressed in [T2], any source 
of disorder will lead to similar results. A recent analysis of a similar model[28j shows that 
the disorder caused by competitive interactions also leads to a resonance effect. 

3 Numerical results 

We have ran the model using different topologies for the neighbourhood network, namely: 
a two-dimensional lattice with von Neumann neighbourhood (and periodic boundary con- 
ditions); a fully-connected network; and a small- world network [29]. In each case, we set 
the amplitude of the forcing A to a small, sub-threshold, value such that, in the absence of 
diversity, the average opinion makes small oscillations around the value m = +1 or m = - 1. 
As the diversity increases, the amplitude of the oscillations first increases and then decreases 
again. 

For the two-dimensional lattice, we plot in figure Q] the spectral amplification factor R 
as a function of diversity a for two different values of the parameter a, which determines 
the relative importance that preference has with respect to the external signal as compared 
to the neighbours influence. Although there are some clear finite-size effects, a well defined 
maximum is clearly observable in both cases. Comparing panels QJa) andQJb), it is apparent 
that for a — 1 the response at the optimal diversity level is lower that for a — 0. The 
reason for this dependency on the parameter a can be easily understood by inspection of 
the dynamic rules of the system (cf. Eq. [3]): this parameter can be seen as a secondary 
source of disorder in the system disturbing the external signal, lowering its effect. Thus, 
the effective level of disorder introduced by the diversity in the population is increased, for 
increasing values of a. This can be confirmed by the fact that in panel (b), the location of 
the maximum response is shifted towards lower values of diversity, measured in terms of a. 
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Figure 1: Spectral amplification factor, R, as a function of the diversity a when the neighbourhood 
network is a two-dimensional regular lattice. Each panel correspond to different values of the 
parameter a, that measures the relative weight of the individual preferences with respect to the 
external signal: (a) a = 0, (b) a = 1. The symbols represent different system sizes: N = 5 2 (o), 
N = 10 2 (□), N = 30 2 (o) and N = 100 2 (A). All the curves show an optimum response for 
an intermediate value of diversity. In both panels, the signal has an intensity A = 5 x 10~ 3 and 
frequency O = 27r/1024. 

Figure [2] shows the results for the small- world network. This is constructed in the usual 
way [29], with a rewiring probability p and average connectivity 2k. The main result here 
is that as p increases (leading to a larger degree of disorder), the resonance peak narrows 
and the optimal response increases. Thus, for slightly disordered networks the systems reacts 
more robustly amplifying the external stimulus for a wider range of diversity values; however, 
R reaches lower values. This effect can be understood if we consider that for small values 
of p, it is more probable for the system to develop stable domains of different opinions, 
that enlarge or shrink depending on the instantaneous signal value. These domains tend to 
coexist, thus decreasing the global response of the system. 

Finally, figure [3] shows the response of the system for a fully-connected neighbourhood 
network. In panel (a), different symbols correspond to different signal amplitudes. We have 
also included in this figure the case of a signal amplitude A — 0.60 which is supra-threshold 
in the zero-diversity case. As seen in this panel, supra-threshold signals are not amplified at 
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Figure 2: For a small- world neighbourhood network, we plot the response R as a function of 
diversity, <r, and the rewiring probability p. It is apparent that, regardless the exact network 
topology, the effect of diversity-induced resonance appears in the system, showing the existence 
of an optimum synchronisation between the external signal and the global dynamics of the 
system. The signal has an intensity A = 5 x 10~ 3 and frequency ft = 27r/1024. The system size 
is N = 10 3 . The initial network is a one-dimensional one with k = 3, i.e. before rewiring, each 
site is connected to its six nearest neighbours. 

all, and the response of the system monotonically decreases with diversity. This trait was 
found in diversity-induced resonance in bistable systems [12], but is also commonly found 
in systems exhibiting stochastic resonance [4] with respect to noise intensity. Also, in the 
same panel, it can be seen that larger (but still sub-threshold) signal amplitudes, lead to 
lower optimal diversity values. In panel (b) of figure [3] we plot the dependence of the system 
response R with respect to the signal frequency f2. As shown, the larger the frequency, the 
lower the response of the system. In this discrete model, this is because slower signals allow 
units with a given bias more akin to adopt its favoured opinion. 

4 Analytical approach 

We now present a mean-field like theory that can explain the observed features. The deriva- 
tion follows the lines of [28j. Since the opinion changes at each time by the modification of 
a single variable, we can write the following exact relation for the ensemble average m(t): 



Nm{t + dt) = Nm(t) + (^(t + dt) - (a(t)\{(i(t)}) 



(6) 



7/02 



Claudio J. Tessone and Raul Toral: 
Diversity-induced resonance in a model for opinion formation 
submitted 



(a) 



10' 



R 



10 



10' 



10, 



1 1 1 1 1 1 1 — 


• > i 1 * 1 * i * 1 * 1 . 


< 


> rj 









(b) 



10 



10 



0.0 0.5 1.0 1.5 2.0 




Figure 3: We compare the theoretical prediction of the system response (see inline text for details) 
with numerical simulations in fully-connected networks. In panel (a) each symbol represents the 
response R, as a function of diversity a for different values of signal intensity: A = 0.01 (•), 
A = 0.05 (□), A = 0.10 (o) and A = 0.60 (A). The full lines correspond to the analytical 
prediction. The other parameters are: N = 10 4 , Q = 27r/1024 and a = 0.1. Panel (b) shows the 
spectral amplification factor as a function of signal frequency, f2 for two different values of the 
signal intensity: A = 0.05 (•) and A = 0.1 (□). The solid lines stand for the theoretical results. 
The other parameters are a = 0.85, a = 0.1 and N = 10 4 . 



where {fi(t)} = (fii{t), . . . , Jijv(£)) denotes the particular realisation of the /i, variables and 
(. . . | . . .) denotes a conditional ensemble average. By identifying dt = 1/N and rearranging 
we get: 

dm(t) 



dt 



-m+(i*(t + dt)\{ii{t)}) 



(7) 



The average term in the right-hand side can be computed using the contribution from the 
biased majority rule, Eq. ©which acts with probability 1— (we define /(£) = Asin(fM)) 
and the contribution of the external forcing, Eq. ([3]) which acts with probability \f(t)\. In 
the spirit of the mean-field approximation we replace in Eq. the average opinion of the 
neighbourhood n(i) by the global average opinion m(t). This yields: 



(m(t + <tt)\{n(t)}) = 



(l-|/(t)|)<sign(m(t)+^)|{MW})+ 
l/^Ksignt/W + a^KM*)}) 



(8) 



8M 



Claudio J. Tessone and Raul Toral: 
Diversity-induced resonance in a model for opinion formation 
submitted 



Both mean values can now be easily evaluated: 

(sign(m(f) + Oi)) = Prob(0; > -m(t)) ~ Prob(6». t < -m(t)) 

= G(m(t)). (9) 

Here, G{6) = 1 — 2G(—0), where G{9) is the cumulative probability function of the distribu- 
tion of preferences g{9). In the same way, the contribution of the external signal to equation 
© is 

( S ign(f(t)+a9i))=G(f(t)/a). (10) 

Adding up those contributions we get a closed evolution equation for the average opinion 
m(t): 

^ = -m + \f(t)\G (/(*)/«) + (1 - |/(t)|) G(m). (11) 

This equation can be written as a relaxational dynamics [30] in a time-dependent potential 

dm(t) _ dV(m,t) 
dt dm ' 1 ' 

It is easy now to see the effect that the diversity has on the dynamics of the global variable. 
For the sake of concreteness, we consider that the preferences follow a Gaussian distribution 
of zero mean and variance a 2 , but similar results hold for other distributions. In this case, 
G{6) = erf(6»/CT V / 2), being erf(x) the error function [3 lj . Consider first the non-forced case, 
fit) = 0. The potential is: 

W = ^- m erf(-^)-,yfe-^. (13 ) 

This potential is bistable when a — 0. As a increases, the two minima of the potential 
get closer to each other and the barrier between them decreases until at the critical value 
a c = \j2j-K the potential becomes monostable, as can be seen in figure [H This shows the 
existence of a phase transition between consensus and non-consensus states. The effect of the 
external field now is easily understood as a periodic lowering and rising of the two potential 
wells. For small a, the barrier separating the two stable points is large and the effect of the 
field is that of making the global variable m(t) oscillate around one of the equilibrium points. 
As a increases, the barrier lowers and it is possible to induce transitions between the two 
stable states. When a increases even further, the potential becomes monostable and again 
the effect of the forcing is that of producing small oscillations, this time around the only 
equilibrium point. This macroscopic mechanism is similar to the one found in |12j . 

We can now integrate numerically equation [TT] to obtain the time evolution of m(t) and 
compute the spectral amplification factor (cf. Eq. [5]) for a given set of parameters. In both 
panels of figure[3l we display with lines this theoretical prediction. As expected from a mean- 
field type theory, there is a good agreement with the numerical results of the fully connected 
network. 
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Figure 4: Effective potential (cf. Eq. [T3|) defining the relaxational dynamics of the average opinion 
for different values of diversity: a = 0, 0.5, 0.7979, and 1 (respectively, from top to bottom as 
found at the vertical line m = 0). 

5 Conclusions 

We have analysed theoretically and numerically a model for opinion formation. The model 
has many points in common with random field Ising models used to study phase transitions 
in statistical mechanics. It incorporates two basic ingredients for the evolution of the opinion 
held by an individual: social pressure and the effect of advertising (modelled as an oscillating 
influence acting over all the individuals). The model also considers that every individual 
has an intrinsic preference for one or the other option. We have shown that an optimal 
synchronisation of the average opinion with respect to the external signal can be achieved 
if the population shows some degree of diversity in the preferred opinions. We have also 
shown that the results are robust against the exact topology of the network used to model 
the neighbourhood of the individuals, and that the results hold for increasingly large system 
sizes. We have given explanations for this resonance both from the point of view of the 
individual responses to the external influence or by looking and the average global variable 
within a mean-field approach. 

From the point of view of the social dynamics, our results imply that an external message 
can propagate better in a society if there is some degree of diversity in the individual pref- 
erences. These results can also be interpreted in the context of population dynamics where 
the external signal stands for a changing environment [32| [33] . The value of the diversity 
parameter specifies to which external condition an individual is best fitted to. In this setting, 
the response is directly related to the average fitting of the population. Within this interpre- 
tation, the results reported in this paper imply that intermediate values of diversity cause a 
better fit of the population to the changing environment. 
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